EPJ manuscript No. 

(will be inserted by the editor) 



OS 

c 
a 
>—» 

(N 



> 
oo 

(N 

o 

On 
On 

B 

C 

o 
o 



13 



Convective and Absolute Instabilities in the Subcritical 
Ginzburg-Landau Equation 

Pere Colet, Daniel WalgraeP, and Maxi San Miguel 

Institute Mediterraneo de Estudios Avanzados, IMEDEA b (CSIC-UIB), Campus Universitat Illes Balears, E-07071 Palma de 
Mallorca, Spain. 

11 December 1998 

Abstract. We study the nature of the instability of the homogeneous steady states of the subcritical 
Ginzburg-Landau equation in the presence of group velocity. The shift of the absolute instability threshold 
of the trivial steady state, induced by the destabilizing cubic nonlinearities, is confirmed by the numerical 
analysis of the evolution of its perturbations. It is also shown that the dynamics of these perturbations 
is such that finite size effects may suppress the transition from convective to absolute instability. Finally, 
we analyze the instability of the subcritical middle branch of steady states, and show, analytically and 
numerically, that this branch may be convectively unstable for sufficiently high values of the group velocity. 

PACS. 47.20. Ky Fluid Dynamics: Nonlinearity - 47.54.+r Fluid Dynamics: Pattern selection; pattern 
formation - 05.40+j Statistical Physics and Thermodynamics: Fluctuation phenomena, random processes, 
and Brownian motion 

1 Introduction tions. Such instabilities lead to the formation of various 

kinds of spatio-temporal patterns jj]]. Well known examples 

Several physico-chemical systems driven out of cquilib- . . . . . . 

are: Kayleigh-Benard instabilities m Newtonian fluids, bi- 

rium present stationnary instabilities of the Turing type, . jtj a, 

nary mixtures, or viscoelastic solutions PU3I, electrohydro- 

or oscillatory instabilities corresponding to Hopf bifurca- . . ...... ..... ,0 „ . 

dynamic instabilities m nematic liquid crystals M , Turing 

a Permanent address: Centre for Non-Linear Phenomena 

instabilities in nonlinear chemical systems 13], convective 

and Complex Systems, Universite Libre de Bruxelles, Campus 

instabilities in Taylor-Couette devices ||, etc. Close to 

Plaine, Blv. du Triomphe B.P 231, B-1050 Brussels, Belgium 
b Electronic address: |http : //www . imedea . uib . es/PhysDept/] such instabilities, the dynamics of the system may usu- 
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ally be reduced to amplitude equations of the Ginzburg- stability threshold, and to analyze the transition from 



Landau type, which describe the evolution of the patterns convective to absolute instability P ]lC| , |ll| , |l5 16 17| , |l8| . [I9f . 

that may appear beyond the bifurcation point 0. However, in the case of subcritical bifurcations, the non- 

. n . , . , , , linearities are destabilizing, which leads to the failure of 

According to the system under consideration, and to 

iii • i ii i T i linear instability criteria. In a qualitative analysis based 

the nature of the instability, these Gmzburg-Landau equa- 

n -inn on the potential character of the real subcritical Ginzburg- 

tions may contain mean now terms induced by group ve- 

, ., . T ,, . r ■ n j n Landau equation, Chomaz |2C(| argued that the transition 

locities. In this case, pattern formation crucially depends n ' v ° 

,, ,. , , , r , n . , , .t, between convective and absolute instability of the trivial 

on the convective or absolute nature of tne instability. 

T , , , r , , , steady state should occur at the point where a front be- 

Let us recall that, when the relerence state is convectively J 1 

, . . , . i.i tween the rest state and the nontrivial steady state is sta- 

unstable, localized perturbations are driven by tne mean J 

n . , , , , , i . , i e tionary in a frame moving with the group velocity. This de- 

now m such a way that they grow m tne moving reference 

r , , , r, , , , . „ ,, fines the nonlinear convective- absolute instability thrcsh- 

trame, but decay at any fixed location. (Jn the contrary, 

, , , . , i . , i.j u old, above which nonlinear global modes are intrinsically 

m the absolute instability regime, localized perturbations ° 

n , , , . iO mi ii • r ii sustained by the dynamics, as discussed by Couairon and 

grow at any hxed location Ihe behavior ol the sys- J J 

. ,., .. , , . i . T Chomaz [El|. This argument relies on the existence of a 

tern is thus qualitatively very dinerent m both regimes. In i^- 11 ° 

,i . , ,,1 . i i ■ ■ j * , unique front between the basic and the bifurcating states, 

the convectively unstable regime, a deterministic system ^ ° ' 

, j i , i j i as it is the case in the subcritical domain where both basic 

cannot develop the expected patterns, except m particular 

, i , ., . , i . . and bifurcating states are linearly stable. In the supercrit- 

experimental set-ups, while m a stochastic system, noise is 

,. i, ^■n ^ 1 • • • i • 1 j ical domain, where the basic state is linearly unstable, or 

spatially amphhed and gives rise to noise-sustamed struc- 

inimFmi a ,i ^111^1 ± 11 in the complex Ginzburg-Landau equation, this front is 

tures ||9|, |10| , [11| . On the contrary, 111 the absolutely unstable * ° mi 

. -,, i i ,1 i 1 not unique any more, and, as commented by van Hecke et 

regime, patterns are intrinsically sustained by the deter- ^ J ' ' J 

. . ,. i . , . , . , ,, , 1 al. |23l, one has to know which nonlinear front solution is 

mimstic dynamics, which provides the relevant selection 

. irrnrr^ii „ , c selected, to determine the nonlinear stability properties of 

and stability criteria 12 13j. Hence, the concepts ol con- J ^ r 

1 * ii i 1 • 1 1 , • i 1 i 1 the basic state, 
vective and absolute instability are essential to understand 

the behavior of nonlinear wave patterns and their stability „,.,,. ,,. , , . . , 

1 J Within this context, our aim m this paper is to con- 

§0- 

tribute to the study of this problem addressing some as- 

The nature of the instability of the trivial steady state pects of it that so far have not been considered. Addi- 

has been studied, cither numerically, analytically and ex- tionally we study stochastic effects. A first aspect con- 

perimentally: In the case of supercritical bifurcations, lin- cerns finite size effects and their influence on the transi- 

ear criteria are appropriate to determine the absolute in- tion from convective to absolute instability for the trivial 
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steady state. Indeed, our numerical analysis of the evo- one proposed by Thual and Fauve p3fl . It could, further- 

lution of perturbations of this state show that it consists more, provide the last building block needed for the undcr- 

in two stages. The first one is devoted to the building of standing of pattern formation in binary fluid convection, 

a front between this state and the bifurcating one. It is as suggested in Q|. 

during the second stage that this front moves outwards In section ^, we recall the dynamical model. In sec- 

or inwards according to the convective or absolute nature tion |^ we discuss the nature of the instability of the triv- 

of the instability. We will show that, although the abso- ial steady state, and present the results of a numerical 

lute instability threshold may effectively be shifted, due analysis of the problem. In section |J, we show, analyti- 

to nonlinear effects, in agreement with p0| ], the first step cally and numerically, that the subcritical middle branch 

of the evolution is sensitive to the size of the system, and of steady states may be convectively unstable, and may 

this may affect the practical determination of the absolute thus be stabilized by mean flow effects in deterministic 

instability threshold, and even suppress it. systems. Finally, conclusions are drawn in section ||. 

A second aspect is the effect of the group velocity on 

2 The Subcritical Scalar Ginzburg-Landau 

the unstable subcritical branch. In the subcritical domain, 

there is a middle branch of steady states, between the Equation 

trivial and the bifurcating ones. In fact, the nature of the For the gake of simplicity wc will consid er, in this paper, 

instability of this branch in the presence of a group ve- systems described by a scalar order paramcte rlikc variable, 

locity has not been considered so far. In the absence of and where the dynamics is given by the real fifth -order 

group velocity, this branch is absolutely unstable. How- Ginzburg _ Landau equat i n, which may be written, in onc- 

ever, the nature of the instability may be modified in sys- dimensional geometries, as §24 : 
terns with group velocity or mean flow effects. We will 

d t A + cd x A = eA + d*A + vA 3 -A 5 + VCxOM), (1) 

effectively show, analytically and numerically, that this 
unstable subcritical branch may be convectively unstable, 
totally or partly, according to the mean flow intensity. Ef- 
fectively, in deterministic systems, unstable states on this 
branch do not necessarily decay in the presence of group 
velocity, while they may remain long lived in stochastic 
systems. This fact may be of practical importance, since 
it provides an alternative way to stabilize the subcritical 

middle branch, which is qualitatively different from the R± = ^(^ ^ V -1 " 2 + 4e) ■ (2) 



For future reference we have added to the equation a 
stochastic term x(cc,i). This models a Gaussian white 
noise of zero mean and variance given by (x(x, t)x{x' , t')) = 
26(x — x')S(t — t'). In the remainder of this section we con- 
sider the deterministic situation with f = 0. 

Bifurcating uniform steady states A(x, t) = R of this 
equation are well known: 
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The linear evolution of the perturbations p± = A— R± growth of perturbations of the steady state. Hence, a reli- 
around these states is then given by: able stability analysis has to include nonlinear terms. As 

discussed by Chomaz and Couairon (2(],|2l| , the nonlinear 
stability analysis of the trivial steady state relies on its 
response to perturbations of finite extent and amplitude. 
Hence, in the case eq. ([!]), without group velocity (c = 0), 
it is sufficient to consider a front solution joining the 
state at x — > — oo to the i?+ state at x — > +oo. 

In the case of the dynamics given by Eq. (Q), this ve- 
locity, Cf may be calculated exactly (^5|, and is found to 
be (cf. fig. |) 



d t p± + cd xP ± = T2R± Vv 2 + 4e p± +d 2 xP± . (3) 
Hence, in the absence of group velocity, the upper branch 

2 

.R+ exists and is stable for — ^- < e, while the middle 

2 

branch i?_ exists and is unstable for — ^- < e < (cf. fig. 
0). 

This picture is, of course, known to change in the pres- 
ence of a finite group velocity c. Let us first recall the 
linear and nonlinear criteria for convective and absolute 
instability for the trivial steady state A = 0. 



Cf = c 1 = —=(-v + 2\/v 2 + 4e) (for < e < ) 

1 ^3 4 4 ' 

3v 2 

3 Linear and Nonlinear Instability of the c f = c* = 2\/e(for— < e) (7) 



Trivial Steady State 



Note that c* is the linear marginal velocity. 
3.1 Analytical Results If the front velocity is negative, which is the case for 



e < — 3w 2 /16, an isolated droplet of the i?+ state embed- 
ded into the state shrinks, and the state is stable. 
On the contrary, if c/ is positive, which is the case for 
e > — 3u 2 /16, R+ droplets grow, and the state is non- 
linearly unstable. The value e = — 3w 2 /16 corresponds to 
ui — e — ck + k 2 (5) the Maxwell construction of phase transitions in which the 

_ trivial and upper branch have equal stability, 
with k = k' + ik". The usual linear instability criterion M 

When c ^ and v — 1, Chomaz pM showed that, in 



The linear evolution around the trivial steady state ^4 = 
is given by 

d t po + cd x p = e p + <9 2 p (4) 
and the corresponding dispersion relation is 



R^=S>£ = (6) 
uk an the unstable domain (e > — 3ir/16), the instability is non- 



and 5R(w(k)) = gives that the trivial steady state is linearly convective (NLC) whenc/ < c, since, in this case, 

convectively unstable for < e < c 2 /4, and absolutely although expanding, a R + droplet is finally advected out 

unstable for e > c 2 /4. of the system. On the contrary, when Cf > c, the insta- 

However, since the nonlinearities of the dynamics are bility is absolute (NLA), since, in this case, i?+ droplets 

destabilizing, the linear terms may possibly not govern the expand in such a way that they finally invade the system. 
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Hence, on generalizing this argument to arbitrary val- noted. As explained before, the boundary conditions for a 

ued of v, one obtains imposing Cf — c in Eq. (R) that the system of size L were taken as follows: A — at x — for 

transition from convective to absolute instability occurs all times and d x A = at x = L. 

& t ; We only discuss here situations where the nonlinear- 

3 , o 2 ,. ,„ i— s ities dominate over mean flow effects, thus where linear 
e a = —{c + —j=vc - it) (fore < V3i>) 



16' V3 

= -c 2 (fore > V3v) (8) 



instability criterion fails. 



4 



(I) A first case corresponds to c = v = 1. In this 

From this result, it appears clearly that, when group ve- c , , , , , 

> j i & f case, the transition from convective to absolute instability 

locity effects dominate over nonlinear ones (c > yfiv), the ghould Qccur ^ ^ = ^ „ Q 21 Thig . g iUustrated by the 

absolute instability threshold remains the linear one. How- , . n n T n n / \ j n_\ 

■' numerical results presented m tig. y. In ng. £5| (a) and (bj, 

ever, when nonlincaritics dominate (v > -4=), the absolute , , , . . .. , ,. r n , , £ 

' v V3" we show the deterministic evolution ol the field A from 

instability threshold decreases, but remains in the e > , • i i... j a n j c n id 

■' ' random initial conditions around A = U and tor e = U.18. 

domain, when v < cv3- It only becomes negative when mi j i c ^ , . , 

J ° I he data connrm the convective nature of the instability 

v > c\/3. This last case is the one originally considered in ,-i • n , , j-uu-u- c c ± 

& J reflectively, we see, m a hrst stage, the building ol a front 

between the trivial state and the bifurcating one, and, in 

a second stage, this front is advected out of the system. 

3.2 Numerical Analysis q q ^ e con trary, for e = 0.23, the instability is absolute, 

™ , , i i , j i as shown in fig. (c) and (d), where the front moves in 

I he above results have been checked through the numcr- s> rj w v n 

■ n . £ , ■, /Ub. Ti7 -n i i the opposite direction, and the bifurcating state invades 

ical integration of the equation (|1|). We will present here ^ ' ° 

some of the data obtained for systems being initially in ^ e s y s ^ em - 

the trivial steady state, and compare them to the predic- The difference between subcritical and supercritical 

tions obtained from the analytical analysis outlined in the behavior is enlightened in fig. | (e) and (f), where the 

preceeding section. To observe a convective instability we neld evolution has been computed with the same parame- 

consider a semi-infinite system with one of the boundaries ters as in % I ( c ) and ( d ) - except that v has been changed 

anchored to the unstable state A(x = 0) = 0. Experi- from +1 to -1 to simulate supercriticality. In this case, the 

mentally, this boundary condition can be achieved using instability should be convective, since the absolute thresh- 

a negative value for the control parameter e for x < 0. old is e = °- 25 > and the results are in agreement with this 

The numerical integrations have been performed using prediction, 
a finite difference method with a spatial step of 5x = The effect of noise in the regime of convective insta- 

0.05 and time step St = 0.001, except where otherwise bility is presented in fig. 0. The field dynamics has been 
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computed for the same values of the parameters as in fig. spite the fact that the perturbation of the zero state at the 

^ (a), but in the presence of noise of different intensities, initial time is much large (see the figure caption). We note 

The noise intensity has been fixed at £ = 10 -6 in fig. that for e < the evolution would still be slower. We first 

H (a) and at £ = 10~ 14 in fig. || (b). In both cases we observe the formation of the front ( initially moving to the 

observe noise sustained structures: Noise is able to sus- right) and much later, when it reaches the upper branch, 

tain finite field amplitudes (positive or negative, accord- invading the whole syste. Hence, when the characteristic 

ing to the +, — symmetry of the system). Weaker noise length needed for the building of the front is larger than 

induces larger healing length for the pattern. Hence, in the system size, the instability is effectively convective, 

the stochastic case, pattern formation is sensitive to sys- although the system should be in the absolute instabil- 

tem size, since the latter has to be larger than the healing ity regime (in the sense of semi-infinite geometries) . For 

length, for the pattern to be able to develop. the parameters chosen in this example and for a length 

L < 2000 one does not observe the decay of the state 

(2) In a second case, we chose c = 0.5, v = 1.5, and 

A = 0. 

this corresponds to the situation presented by Chomaz 

iF^j , , . a , , . / /7T \ i It is worthnoting that the observed finite size effects 

p(J|| , where nonlinear effects dominate (v > cyo.) and ° 

, r . . , i , . , i ., confirm and complement the analysis made by Chomaz 

where the transition from convective to absolute mstabil- 1 J J 



., . , ... n , . . n m n and Couairon 1 26 of fully nonlinear solutions of Ginzburg- 

lty occurs m the subcritical domain since e a ~ — U.21. lor < — ■ ° 

■ i . , i , i , t , i . . Landau equations in finite domains. In case (1), for e = 

e a < e, the instability is absolute, but the dynamics is w 

,., . , j.rr -r Tin 0.23, nonlinear global (NLG) modes exist, even in finite 

qualitatively different if e is positive or negative. When ' ° v ' ' 

_ r, t, ji- j domains. However, since the basic state, A = is linearly 

e > 0, both linear and cubic terms are destabilizing, and J 

, ., ,. , ,. rr , , , . . , i absolutely stable, NLG modes only develop if the initial 

the building and propagation of fronts between trivial and J 

, .„ . , r , , r „ i condition is sufficiently large for the transients to reach an 

bifurcating states is much taster than for e < 0, when the 

n . , , , .,. . j , , , . j , order one amplitude in the finite domain. Since the am- 

hnear term is stabilizing, and the cubic one is destabihz- r 

TIT1 , , , . , , , plification factor increases exponentially with L, the mini- 

mg. When the dynamics becomes very slow the time and 

, , , i- rr^r mum amplitude of initial perturbations able to trigger the 

system size needed to see the lormation ol a front from 00 

• • . i u ,• i. i ,i NLG mode decreases exponentially with L |2ql. As a re- 

an initial perturbation become very large, so that even m 1 - «— v 

, , , . . . , suit, the development of NLG modes is almost insensitive, 

the absolutely unstable regime one might not observe the 

, A „ . n ., ,. e n ., in most practical situations, to system size, 

decay ol the state A = in finite times for a finite system. 1 ' J 

This effect is illustrated in fig. |B| which corresponds to the On the contrary, in case (2), the basic state is linearly 

absolutely unstable regime. Note the significant increase stable, absolutely and convectively, and the minimum am- 

of the times scales in comparison with figs. || and |], de- plitude of initial perturbations able to trigger NLG modes 
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in finite boxes decreases linearly with L. It is why, in the the case for the evolution of the perturbations of the trivial 
conditions of our numerical analysis, no global mode is steady state since the bifurcation is subcritical, but it may 
obtained for L < 2000. also be the case for the perturbations around the middle 

steady state branch, whose evolution is given by 

4 Stability Analysis of the Bifurcating States dtp _ + cQxp _ = +2 j£ VtfT^p. + d 2 p _ 



4.1 Analytical Results - R - O - 5vV +4e) p 2 _ - (Av - 5vV +4e) p 3 _ 

- 5R-pt - pi (12) 

The linear evolution around the upper branch steady states 

R + and middle branch steady states i?_ is given by eq. (||). The quadratic nonlinearity is destabilizing for e < e L = 

The upper states R + are linearly stable for all e > -0.21 v 2 . In such cases, one has to perform a nonlinear 

On the other hand, the usual linear instability criterion analysis of the dynamics to determine the convective or 

shows that the R- steady states are convectively unsta- absolute nature of the instability. 



ble for d>R 2 _^/v 2 + 4e < c 2 , and absolutely unstable for In the regime where the nonlinearities of the evolution 

&R 2 _\Jv 2 + 4e > c 2 . In other words, these states are abso- equation (|l|) are stabilizing, i.e. for e L = -0.21 v 2 < e < 

lutely unstable in the range 0, the results of the linear analysis may be assumed to be 

I c 2 c 2 valid. Hence, we may safely rely on these results above 

— — (v 2 + — +!)v v 2 — c?) < e < — g( l ' 2 + "^ — u v v 2 — c 2 ) 

the metastability point, i.e. for cm — — 3/16u 2 < e < 0. 

(9) 

Below the metastability point, i.e. for — 0.25w 2 < e < 

and convectively unstable in the windows defined by 

— 3/16v 2 , one has to perform a nonlinear analysis, which, 



y 2 1 2 c 2 

4~ < 6 ^ ~ g( w ~~2 v v v ~~ c ) (10) in this case, relies on the evolution of fronts between mid- 



an j die branch states and the trivial steady state. We do not 

_ }_ f y 2 _|_ ^_ _ yy J v 2 _ c 2) < e < o (11) perform this analysis here since it would only affect quanti- 

8 2 

tatively but not qualitatively the results presented above. 

Hence, when v 2 < c 2 , these steady states are always 
linearly convectively unstable. Still, when v 2 > c 2 , there 

4.2 Dynamics of the Subcritical Unstable Branch 

is a range of linear absolute instability in the middle of 

their domain of existence, and a range of linear convective We have numerically confirmed the convective nature of 
instability close to the points where these states disappear, the instability of the subcritical middle branch. The nu- 
This is shown in fig. |^. merical integration has been performed as indicated in 



Nevertheless the linear stability criteria may fail in the subsection 3.2. Also, as indicated in that subsection, to 



presence of destabilizing nonlinearities. This is not only observe a convective instability we consider a semi-infinite 
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system with one of the boundaries anchored to the unsta- is a range of parameters in which it is only convectively 

ble state. Here we have to take A(x = 0) = i?_ corre- unstable. 

sponding to the field amplitude of the subcritical middle The effect of noise in the convectively unstable regime 

branch. Experimentally, this boundary condition can not of the trivial state was to sustain a structure continuously 

be achieved as easily as before because there is no value excited by noise. In the case of the middle branch. 

of the control parameter e for which A(x) = R- is an ho- and when this is convectively unstable, noise forces the 

mogeneous steady stable state. However depending on the system to relax randomly to either of the two coexisting 

system it can be imposed in different ways. In an optical stable branches, as shown in fig. |^. Still, if noise is weak 

system, for example, the left boundary condition could be in comparison with the strength needed to see its effect in 

achieved injecting an external field at x — with the apro- a finite system, one would observe the middle branch as 



priate amplitude. Finally, as in subsection 3.2, the right effectively stable, 
boundary condition is taken as d x A — at x = L 

We computed the evolution from an initial steady state ^ Conclusions 

withi? 2 . = 0.1 on the middle branch, which corresponds to in t ^[ s pa per, we considered systems described by the 

v = 1 and e = -0.09. Wc then study the system dynamics subcritical Ginzburg-Landau equation, and analyzed some 

for different values of the group velocity c. According to problems related with the effect of group velocities on the 

the previous discussion, for stability of its steady states. In the case of the trivial 

,. s _ „ _ . . . . . . , , , . steady state, it is known that the transition between con- 

(1) c < 0.8, the state K_ should be absolutely unstable 

vective and absolute linear instability regimes is shifted 

(2) c > 0.8, the state i?_ should be convectively un- 

by the effect of destabilizing nonlincaritics, and the cor- 

stable. 

responding nonlinear absolute instability threshold may 

In fig. (a), we present, for c = I, the results obtained easily be computed for semi-infinite systems f2(i| , ^l|| . Our 

for the deterministic evolution of an initial perturbation of numerical study of the evolution of perturbations from the 

the state They show that the instability is effectively trivial steady state in finite systems shows that, in a first 

convective. On lowering the group velocity from c = I step, a front is built between this state and the bifurcat- 

to c = 0.55, the nature of the instability changes from ing one, which corresponds to the upper branch of steady 

convective to absolute, as expected, and shown in fig. [7] states. Then, according to the intensity of the group ve- 

(b). These results confirm that the middle branch, which locity, the front moves outwards or inwards, which corre- 

is always stable for c = 0, may be stabilized by mean flow sponds to convective or absolute instability, respectively, 

effects in deterministic systems, in the sense that there When the characteristic length needed for the building of 
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the front is shorter than the system size, the nature of the 6 Acknowledgements 

instability is in agreement with the theoretical predictions 



made for semi-infinite systems. However, our numerical re- 
sults show that, if the characteristic building length of the 
front is larger than the system size, one will never see in- 
ward motion of the front, and, in this case, even above the 
absolute instability threshold, the instability is effectively 
convective. 

We also studied the instability of the subcritical mid- 
dle branch of steady states, a problem that had not been 
addressed up to now. It may be shown, already at the level 
of a linear analysis, that this branch, which is absolutely 
unstable without group velocity, may entirely become con- 
vectively unstable in the presence of group velocities larger 
than some well-defined critical value. This result has been 
confirmed by the numerical analysis of the evolution of 
perturbations of steady states on this branch. The stabi- 
lization of such steady states has effectively be obtained, in 
deterministic systems, for group velocities in the predicted 
range. In stochastic systems, however, these steady states 
relax to one of the stable branches, as expected. Neverthe- 
less, for this relaxation to occur, either noise strength or 
system size have to be large enough. This effect may be of 
practical importance, for example, in binary fluid convec- 
tion, where, besides the fact that the role of subcriticality 
is not clearly understood yet |^7j , the presence of natural 
or forced mean flows, or group velocities, could effectively 
stabilize otherwise unstable branches of steady states. 
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Fig. 3. Deterministic evolution of perturbations of the trivial steady state. The column on the left shows the spatio-temporal 
evolution of the field A and the column on the right the value of the field at different times. The initial condition for each grid 
point is Ai(t — 0) = f\iji\, where r) is a Gaussian random number of zero mean and variance 1 and / = 1CT 4 . In (a) and (b) 
we consider the subcritical case in a convectively unstable regime, with e = 0.18, v = 1, c = 1. In (c) and (d) we consider the 
subcritical case in an absolutely unstable regime, with e = 0.23, v = 1, c = 1. In (e) and (f) we consider the supercritical case, 
with e = 0.23, v = — 1, c = 1. Dashed lines in (b) and (f) correspond to early times when the front is being formed and it 
effectively moves to the left. Continuous lines show the front moving to the right in (b) and (f) and moving to the left in (d) 
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Fig. 4. Spatio-temporal evolution of the field A with noise 
in a convectively unstable regime for e = 0.18, v = 1, c = 1. 
The initial condition is A(x, t = 0) = 0. (a) noise intensity 
£ = 10" 6 , (b) noise intensity £ = 10" 14 . 



Fig. 5. Spatio-temporal evolution of the field A without noise 
in an absolutely unstable regime with e = 0, v = 1.5, c = 0.5. 
The initial condition is as in Fig. ^| but with / = 10 -2 . In this 
case we have taken 8x = 1, St = 0.2. 

2 



convectively unstable 




-0.25 -0.21 e/v 2 

Fig. 6. Domains of linear convective and absolute instability 
of the middle branch uniform steady state in the ((c/w) 2 , e/v 2 ) 
plane. The linear stability analysis is not valid in the hatched 
domain where the nonlinearities are destabilizing. 
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Fig. 7. Deterministic evolution of a perturbation around the 
middle branch state R- in a convectively unstable regime for 
e — —0.09, v — 1. The initial condition for each grid point is 
Ai(t = 0) = R- +/|»7i|, where r\ is a Gaussian random number 
of zero mean and variance 1 and i?_ = VO.l. (a) convectively 
unstable regime with c = 1 and / = 10~ 2 . (b) absolutely un- 
stable regime with c = 0.55 and / = 10~ 6 . 



